Efficient calculation of the Greeks for exponential Lévy processes: an application of measure valued differentiation

Monte Carlo simulation methods have become more and more important in the financial sector in the past years. In this paper, we introduce a new simulation method for the estimation of the derivatives of prices of financial contracts with respect to (w.r.t.) certain distributional parameters called the ‘Greeks’. In particular, we assume that the underlying financial process is a Lévy-type process in discrete time. Our method is based on the Measure-Valued Differentiation (MVD) approach, which allows representation of derivatives as differences of two processes, called the phantoms. We discuss the applicability of MVD for different types of option pay-offs in combination with different types of models of the underlying and provide a framework for the applicability of MVD for path-dependent pay-off functions, as Lookback Options or Asian Options

Behavioral pricing of energy swing options by stochastic bilevel optimization

Holders of energy swing options are free to specify the amounts of energy to be delivered on short notice, paying a fixed price per unit delivered. Due to the complexity of potential demand patterns, risk elimination by replication of these contract at energy exchange markets is not possible. As a consequence, when selling delivery contracts, the energy producer has to explicitly consider the risk emanating from fluctuations in supply cost. The impact of these risk factors can be mitigated by the contract seller, who is an energy producer, to a certain extent: Supply cost fluctuations can be absorbed by the own generation portfolio whereas demand uncertainties can be influenced by the choice of the strike price, implicitly changing the buyer's behavior. Considering this, the determination of the optimal strike price can be formulated a a stochastic bilevel problem where the optimal decision of upper level player (price setting and production) depends on the optimal decision of a lower level player (demand depending on the price). We present a solution algorithm tailored to the resulting specific stochastic bilevel problem. We illstrate the effects of the behavioral pricing approach by studying behavioral price setting for natural gas swing options, highlighting in particulr the influence of the seller's production and contract portfolio as well as of the market liquidity on optimal exercise prices

Guaranteed Bounds for General Nondiscrete Multistage Risk-Averse Stochastic Optimization Programs

In general, multistage stochastic optimization problems are formulated on the basis of continuous distributions describing the uncertainty. Such “infinite” problems are practically impossible to solve as they are formulated, and finite tree approximations of the underlying stochastic processes are used as proxies. In this paper, we demonstrate how one can find guaranteed bounds, i.e., finite tree models, for which the optimal values give upper and lower bounds for the optimal value of the original infinite problem. Typically, there is a gap between the two bounds. However, this gap can be made arbitrarily small by making the approximating trees bushier. We consider approximations in the first-order stochastic sense, in the convex-order sense, and based on subgradient approximations. Their use is shown in a multistage risk-averse production problem

Optimal XL-insurance under Wasserstein-type ambiguity

We study the problem of optimal insurance contract design for risk management under a budget constraint. The contract holder takes into consideration that the loss distribution is not entirely known and therefore faces an ambiguity problem. For a given set of models, we formulate a minimax optimization problem of finding an optimal insurance contract that minimizes the distortion risk functional of the retained loss with premium limitation. We demonstrate that under the average value-at-risk measure, the entrance-excess of loss contracts are optimal under ambiguity, and we solve the distributionally robust optimal contract-design problem. It is assumed that the insurance premium is calculated according to a given baseline loss distribution and that the ambiguity set of possible distributions forms a neighborhood of the baseline distribution. To this end, we introduce a contorted Wasserstein distance. This distance is finer in the tails of the distributions compared to the usual Wasserstein distance

Assessing current and future impacts of climate-related extreme events. The case of Bangladesh

Extreme events and options for managing these risks are receiving increasing attention in research and policy. In order to cost these extremes, a standard approach is to use Integrated Assessment Models with global or regional resolution and represent risk using add-on damage functions that are based on observed impacts and contingent on gradual temperature increase. Such assessments generally find that economic development and population growth are likely to be the major drivers of natural disaster risk in the future; yet, little is said about changes in vulnerability, generally considered a key component of risk. As well, risk is represented by an estimate of average observed impacts using the statistical expectation. Explicitly accounting for vulnerability and using a fuller risk-analytical framework embedded in a simpler economic model, we study the case of Bangladesh, the most flood prone country in the world, in order to critically examine the contribution of all drivers to risk. Specifically, we assess projected changes in riverine flood risk in Bangladesh up to the year 2050 and attempt to quantitatively assess the relative importance of climate change versus socio-economic change in current and future disaster risk. We find that, while flood frequency and intensity, based on regional climate downscaling, are expected to increase, vulnerability, based on observed behaviour in real events over the last 30 years, can be expected to decrease. Also, changes in vulnerability and hazard are roughly of similar magnitudes, while uncertainties are large. Overall, we interpret our findings to corroborate the need for taking a more risk-based approach when assessing extreme events impacts and adaptation - cognizant of the large associated uncertainties and methodological challenges -

Modeling dependent credit rating transitions: a comparison of coupling schemes and empirical evidence

Three coupling schemes for generating dependent credit rating transitions are compared and empirically tested. Their distributions, the corresponding variances and default correlations are characterized. Using Standard and Poor's data for OECD countries, parameters of the models are estimated by the maximum likelihood method and MATLAB optimization software. Two pools of debtors are considered: with 5 and with 12 industry sectors. They are classified into two non-default credit classes. First portfolio mimics the Dow Jones iTraxx EUR market index. The default correlations evaluated for 12 industry sectors are confronted with their counterparts known for the US economy

Incorporating model uncertainty into optimal insurance contract design

In stochastic optimization models, the optimal solution heavily depends on the selected probability model for the scenarios. However, the scenario models are typically chosen on the basis of statistical estimates and are therefore subject to model error. We demonstrate here how the model uncertainty can be incorporated into the decision making process. We use a nonparametric approach for quantifying the model uncertainty and a minimax setup to find model-robust solutions. The method is illustrated by a risk management problem involving the optimal design of an insurance contract

Numerical Modeling of Dependent Credit Rating Transitions with Asynchronously Moving Industries

Two models of dependent credit rating migrations governed by industry-specific Markovian matrices, are considered. Caused by macroeconomic factors, positive and negative unobserved tendencies, encoded as values "1" or "0" of the corresponding variables, modify the transition probabilities and render the evolutions dependent. They are neither synchronized across industry sectors, nor over credit classes: an upswing in some of them can coexist with a decline of the rest. The models are tested on Standard and Poor's data. MATLAB optimization software and maximum likelihood estimators are used. Obtained distributions of the hidden variables demonstrate that the considered industries migrate asynchronously trough credit classes. Since downgrading probabilities are less affected by the unobserved tendencies, estimated by Monte-Carlo simulations distributions of defaults, exhibit lighter, than for the known coupling models, tails for schemes with asynchronously moving industries. Moreover, the lightest tails were obtained in the case of industry-specific transition matrices

Critical growth of cerebral tissue in organoids: theory and experiments

We develop a Fokker-Planck theory of tissue growth with three types of cells (symmetrically dividing, asymmetrically dividing and non-dividing) as main agents to study the growth dynamics of human cerebral organoids. Fitting the theory to lineage tracing data obtained in next generation sequencing experiments, we show that the growth of cerebral organoids is a critical process. We derive analytical expressions describing the time evolution of clonal lineage sizes and show how power-law distributions arise in the limit of long times due to the vanishing of a characteristic growth scale